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ABSTRACT 

Accurate and complete rotational, rotational-vibrational and rotational-vibrational- 
electronic line lists are calculated for sodium hydride: both the NaH and NaD iso- 
topologues are considered. These line lists cover all ro-vibrational states of the ground 
(X ^11+) and first excited {A ^I]+) electronic states. The calculations use available 
spectroscopically-determined potential energy curves and new high-quality, ab initio 
dipole moment curves. Partition functions for both isotopologues are calculated and 
the effect of quasibound states is considered. The resulting line lists are suitable for 
temperatures up to about 7000 K and are designed for studies of exoplanet atmo¬ 
spheres, brown dwarfs and cool stars. In particular, the NaH A — X band is found 
to show a broad absorption feature at about 385 nm which should provide a signa¬ 
ture for the molecule. All partition functions, lines and transitions are available as 
Supplementary Information to this article and at www.exomol.com, 

molecular data; opacity; astronomical data bases: miscellaneous; planets and satellites: atmospheres; stars: low-mass 


1 INTRODUCTION 

Although long discussed ( Kirby fc DalgarriJ 1978ll . the existence of sodium hydride, NaH, is yet to be confirmed in any 
astronomical bodies. Notably few metal h ydrides have been detected in the inters t ellar medium and s earches for NaH in both 
dense clouds ( Plambeck fc EricksonlllQsT and diffuse clouds | Snow fc Smith|[l977 ; Czarnv et al.lll9^ ) have so far only yielded 
upper limits. Similarly an a ttempt to detect the long-wavelength signature for NaH in the atmosphere of Jupiter and Saturn 
i Weisstein fc Serabvn 19961) also proved negative. _ 


While spectroscopic data for lon g-wavelength studies of NaH has been the su bject of the laborat ory studies llLeopold et al. 


19871: lOkabavashi &; Tanimoto 


and are compiled in the CDMS database ( MulleretaU|20& 


the situation is less clear 

at infrared and visible wavelengths. A recent work on M-dwarf models by iRaipurohit et al.l ll2013t) identified NaH as one of 
only three molecules for which the necessary spectroscopic data is missing. This is despite the fact that the A — X band of 
NaH should give an observable feature at visible wavelengths. It is this situation we address here. We note th at the next 
lowes t-lying singlet electronic excited states of NaH are the C i Walii et al. 20151) and the shallow B ^H state ( Yang et al. 


2004 ) whose yet to be observed spectra lie in the ultraviolet, probably at around 3000 A. 


The large dipole moments of NaH, as well as supposedly ma king it more amenable to astronomical detection, has also 
excited the interest of scientists working on ultra-cold molecules ( Juarros et al. 20061 : Avmar et al. 20091 : Su et al. 20ld) . In 
particular, it has been sug gested that these dip oles should enhance the prospects of molecule formation by radiative association 


at very low temperatures (lJuarros et al.l l2i 


The ExoMol project aims to provide line lists of spectroscopic transitions for key molecular species which are likely to be 
im portant in the atmospheres of extrasolar planets and cool stars; its aims, scope and methodology have been summarised 
by jTennvson fc Yurchen^ ( 20121 ). In this paper, ro-vibrational and rovibronic transition lists are computed for the only two 
stable isotopes of sodium hydride, ^®NaH and ^®NaD. These line lists consider transitions within and between the two lowest 
electronic states of NaH, X and A ^E"*". The resulting line lists are comprehensive and should be valid up to the 7000 K 
temperature range. 
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2 Rivlin et al 



Figure 1. Empirical NaH and 


state potential energy curves 


due to IWalii et alJ l|2Q15h . 


2 METHOD 


As w i th other diatomic syst ems studied as part of the ExoMol project ( Yadin et al hold : iBarton et al.ll20 13l.l2014l : lYorke et al. 


2014 : iPatrascu et al.N2015l l. we treated NaH using both ab initio and experimental data. IWalii et al.l (120151 1 used a direct 


potential fit analysis to generate accurate potential energy curves for the ground and first singlet excited states of NaH 
{X and ‘Avoided Cro ssing’ A ^E^). These curves (also shown In Flg.[T| were then used as input to the program LEVEL 
version 8.2 (|Le R,ovll2007l l to produce accurate nuclear motion energy levels, wavefu nctions and hence t ransition intensities. 


A full discussion of how the potential energy curves were derived can be found in Walii et al. 1 2015l l. which also gives a 


compre hensive survey of p revious laboratory work on the spectrum of NaH and NaD. We note that the experimental data 
used bv lWalii et al.l ( 2015l l to determine their curves is insufficient to char acterise the last 10 00 cm“^ of the excited A ^E^ 
below dissociation, where the A state undergoes another avoided crossingj|Ajmaret_^ 200(^ 1. The dissocation limit of this 
state was constrained by the atomic limit neglecting spin-orbit coupling I Walii et’aH^OlS l and ro-vibrational states computed 
in this region must be considered approximate. 


2.1 Dipole Moments 


Ab initio dipole moments were computed on a dense, uniformly-spaced grid of 220 points from r=2.40 to r=13.40 ao in steps of 
0.05 ao. Dipoles were computed as expectation value of a multi-reference configuration interaction (MRCI) calculation using a 
full valence complete active space model, which distributes the two outer electrons among 5 orbitals, and the aug-cc-pwCV5Z 
basis set. Orbitals were produced using a state-averaged complete active space self consistent field (CASSCF) which considered 
both the X and A states. The core 2s2p electrons were correlated at the MRCI stage. A single run takes about 4 minutes. 

Figure [2] shows our two diagonal dipole moment curves and the off-diagonal A-X transition moment. The diagonal dipoles 
have a somewhat unusual shape which is undoubted associated with the changing character of the N aH wavefunction as a 
function of geometry. The dipole moment of ground state NaH is has been determined experimentally I Dagdigian 1979I L but 
not particularly accurately. Conversely there have been quite a numbe r of theoretical stud ies on the problem. The most recent 
Is a comprehensive calculation of all three moment considered here bv lAvmar et al.l (|200S ), who als o prov ide a comprehensive 


survey of the previous literature. Our results are in excellent agreement with those lAvmar et al 
various turning points in the dipole curves which are a feature of our calculations. 


( 2 ^), 


who also find the 


Our dipole moment points, which are provided in the supplementary material, were input directly Into LEVEL. 


3 RESULTS 

3.1 Partition Function 


Partit ion functions for NaH and NaD were calculated by summing over all energy levels calculated. ISzldarovszkv fc Csaszar 
1 2 OI 5 I 1 have recently shown that inclusion of quasibound (resonance or pre-dissociative) states can influence the partition 
function sums at higher temperatures. Here we present partition functions for both NaH and NaD obtained by summing both 
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Figure 2. Ab initio NaH dipole moment curves as a function of internuclear separation. 


over all the truly bound states of the system and also by including quasibound states determined using LEVEL. Table [T] 
compares these two methods and suggests that below 2000 K the influence of quasibound states on the partition sum can safely 
be neglected. Above this temperature the quasibound states have a small but increasing influence. However for consistency 
with our previous work and because we do not consider transitions to quasibound states, we ignore the effect of quasibound 
states below. For high accuracy work at higher temperatures this approx imation may need to be re-visited. 

Table [T] also c ompares our NaH partition functions with those of ISauval fc TatumI 111984 1 and the CDMS database 
I Muller et al. 2005h . Our results cover a larger temperature range: CDMS values only go up to 300 K for NaH, and Sauval 
& Tatum only claim their values are accurate for 1000 - 9000 K. The latter values were multiplied by the factor 8 in order 
to explicitly account for the nuclear spin degeneracy in accord with the HITRAN ( Fischer et al. 2003l l convention, used by 
ExoMol; the CDMS partition function had to be multiplied by 2. According with th e convention of CDMS, the common 
factors in partition functions are usually divided to keep them small ( Muller et al. 2005h . After these adjustments results give 
good agreement with Sauval & Tatum and CDMS in the appropriate temperature ranges, although there is a slight divergence 
from Sauval & Tatum’s values at higher temperature ranges. 

A polynomial fit of this work’s partition functions is given in in Table [5] The values in Table [5] are coefficients in the 
following equation: 


logio Q{T) = ^ an(logio T)", 

n=0 

which was based on a series expansion used by Vidler fc Tennyson ( 2000ll . 


( 1 ) 


3.2 Line lists 

Line lists were calculated for the two main isotopologues: NaH and NaD. All transitions satisfying the selection rule A J = ±1 
between all ro-vibrational states in the electronic ground state and and first electronically-excited state were considered. A 
summary of each line list is given in Table [S] Although every possible transition was computed, some were very weak and are 
not included in the final line list. 


4 RESULTS 


The line lists contain hundreds of thousands of transiti ons. They are separat ed into separate energy levels and transitions 
files. This is done using the standard ExoMol format ( Tennyson et al. l2013l i where the states file gives all energy levels 


and associated quantum numbers, and the transitions file gives Einstein A coefficients and the numbers of the associated 
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Table 1. Partition function of NaH compared with other sources 


T(K) 

This work 
(with quasibound) 

This work 
(without quasibound) 

Sauval Tatum f19841 

CDMS 

37.5 

38 

46.50 

46.50 

51.00 

51.67 

45.93 

75 

89.10 

89.10 

97.98 

89.10 

150 

175.73 

175.73 

182.99 

175.72 

225 

262.83 

262.83 

266.93 

262.62 

300 

351.39 

351.39 

353.94 

349.79 

500 

607.74 

607.74 

610.87 


1000 

1479.85 

1479.85 

1460.52 


1500 

2764.91 

2764.89 

2668.35 


2000 

4521.88 

4520.67 

4303.28 


2500 

6811.55 

6799.15 

6437.62 


3000 

9692.63 

9633.66 

9149.42 


3500 

13205.11 

13024.25 

12523.04 


4000 

17363.28 

16941.48 

16649.49 


4500 

22161.35 

21341.29 

21626.79 


5000 

27583.33 

26178.70 

27560.23 


5500 

33610.90 

31415.33 

34562.72 


6000 

40227.09 

37021.39 

42755.06 



Table 2. Fitting parameters for the partition function of NaH/NaD. Fits are valid for temperatures up to 9000 K within 1 %. 



NaH 

(with quasibound) 

NaH 

(without quasibound) 

NaD 

(with quasibound) 

NaD 

(without quasibound) 

ao 

5.20312719799 

-0.19230735479 

-1.61166291473 

-0.19230735479 

ai 

-19.01261024630 

3.70859873334 

9.53281763105 

3.70859873334 

a2 

34.83255367250 

-3.07887332529 

-13.18395704240 

-3.07887332529 

0.3 

-34.37463416770 

0.71923866264 

10.40662225870 

0.71923866264 

a4 

20.73051960960 

1.11929540933 

-4.49578453850 

1.11929540933 

05 

-7.76051252085 

-0.99902883569 

1.01674402973 

-0.99902883569 

aQ 

1.75196812228 

0.34706733327 

-0.09068752170 

0.34706733327 

07 

-0.21734537482 

-0.05612884704 

-0.00355980420 

-0.05612884704 

as 

0.01134093577 

0.00348875707 

0.00081789072 

0.00348875707 


Table 3. Summary of the computed NaH and NaD line lists. 



NaH X state 

NaH A state 

NaD X state 

NaD A state 

Maximum v 

21 

32 

30 

34 

Maximum J 

81 

123 

113 

171 


Table 4. Extract from the states file for NaH. The zero of energy is taken to be the energy of the lowest energy level. The files contain 
3 339 levels for NaH and 5 960 levels for NaD. 


n 

E 

g 

J 

State 

V 

1 

0.000000 

8 

0.0 

X 

0 

2 

1133.102453 

8 

0.0 

X 

1 

3 

2228.214399 

8 

0.0 

X 

2 

4 

3285.975269 

8 

0.0 

X 

3 

5 

4306.940010 

8 

0.0 

X 

4 

6 

5291.562054 

8 

0.0 

X 

5 

7 

6240.146564 

8 

0.0 

X 

6 

8 

7152.840057 

8 

0.0 

X 

7 

9 

8029.619734 

8 

0.0 

X 

8 

10 

8870.226718 

8 

0.0 

X 

9 

11 

9674.057941 

8 

0.0 

X 

10 

12 

10440.079797 

8 

0.0 

X 

11 


n: State counting number. 

E\ State energy in cm“^. 
g: State degeneracy. 

J: State rotational quantum number. 
v: State vibrational quantum number. 
State: Electronic state label. 
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Table 5. Extract from the .trans file for NaH. Full tables are available from the Exomol web site www.exomol.com or from the VizieR 
service from the web site of the Strasbourg Data Center cds.u-strasbg.fr The files contain 79 898 transitions for NaH and 167 224 
transitions for NaD. 


/ 

i 

^if 

342 

322 

2.1771E+01 

3190 

3198 

5.9304E+00 

1168 

1200 

6.8873E-01 

3329 

3327 

1.1725E+02 

477 

495 

5.0419E-01 

3225 

3217 

8.4415E+01 

2887 

2870 

3.5561E+01 

929 

938 

1.0760E-01 

759 

744 

1.1033E+01 

3286 

3280 

9.6810E+01 


/: final state number. 

V. Initial state number. 

Aij: Einstein A coefficient in 


Table 6. Lifetimes (i n ns) of selected ro-vibrational states within the A electronic state: comparison of experiment IPagdigianl 

(Il976h ; iNedeled lll983li and theory (this work). 


v' 

J' 

Ref. 

Exp. 

Calc. 

3 

8 

Dagdiffian (1976) 

24.0(3) 

25.71 

4 

11 

Dagdigian (1976) 

28.3(3) 

25.69 

5 

16 

Dagdigian (’1976) 

27.1(3) 

26.44 

3 

4 

Nedelec Q9831 

27 

25.37 

4 

6 

Nedelec Q9831 

28.5 

25.16 

7 

6 

Nedelec fl983') 

26.5 

25.40 

9 

6 

Nedelec {'19831 

26 

29.71 

11 

6 

Nedelec ('1983') 

25.5 

33.10 

14 

7 

Nedelec ("19831 

27 

36.49 

18 

11 

Nedelec ("19831 

27.5 

44.01 

21 

6 

Nedelec (19831 

35 

48.67 


states. Extracts from the start of the NaH states and transitions files are given in Tables |4] and [5] The full line list can be 
downloaded from the Exomol web site www. exomol. com or from the VizieR service from the web site of the Strasbourg Data 
Center cds .u-strasbg. f r, A small computer program which uses these files to compute spectra is given in the supplementary 


material to Yorke et al 


The direct potential fit Walii et al. 1 20151) essentially reproduced all experimental data within its stated uncertainty oi 

a ve ry modest a moun t of data cleaning had b een performe d. _ 

Dagdigj^ 1 197fil ). Baltavan et al. 1 197(ih and Nedelec 1 198,11 ) measured radiative life times for a number of A state 


vibrational levels. These measurements suggest that the lifetime actually has little dependence on the actual level involved 
with a sudden increase for v = 21. A compariso n with ou r pred ictions for some of these states is listed in Table[6] in general the 
agreement is very good. However, in contrast to Nedelec 1 19831) . we observe the increase of the life time of the A state levels at 
lower vibrational excitations than v = 21. The increase is associated with the X state dissociation threshold of 15797.4 cm“^, 
which truncates the contributions from the bound states of the ground electronic state. Thus the lower experimental values 
of the life times can be attributed to the transitions to the X quasibound states, which are not included in our line lists. 

Figure [3] gives an overview of the entire spectral range covering the X and A systems, from infrared to visible shown as 
an absorption spectrum generated for T = 2000 K. Of particular astronomical interest is likely to be the A — X band in the 
380 nm region. The A — X band does not exhibit sharp features which can be attributed to the shallow A potential energy 
curve, see Fig. [T] This reduces the chances of NaH to be detected at visible wavelengths, especially at high temperature, as 
also illustrated by Fig. |4j where the A — X absorption band is shown for a set of temperatures from 1000 to 4000 K. 

Figure [5] shows the pure rotational band of NaH at T = 298 K compared to the spectra from the CDMS database. The 
latter was generated by averaging over the the hyperfine lines, which are not resolved in o ur line lists. The ag reement between 
intensities can be attributed to the equilibrium X-state dipole moment value of 6.7 D of I Sachs et al.l (119751 ) used by CDMS. 
Our value is 6.43 D at re = 1.887. 
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wavelength, |,im 

0.2 1.0 0.67 0.5 0.4 0.33 



0 5000 10000 15000 20000 25000 30000 

wavenumber, cm ' 


Figure 3. Overview of the X-X, A—A, and A-X absorption spectrum of ^®NaH at 2000 K. The theoretical spectrum is convolved with 
a Gaussian profile with a half-width half-maximum (HWHM) of 0.1 cm“^. 



Wavelength, A 


Figure 4. Computed A-X absorption spectrum of ^^NaH at different temperatures. A Doppler profile was used to convolve the theoretical 
spectrum. Wavenumber bin of 50 cm“^ was chosen to mimic resolving power of about 500 to 1000. 


5 CONCLUSIONS 


We present comprehensive line lists for NaH and NaD. These are based on the direct solution of the nuclear motion Schrodinger 
equation us i ng po tential energy curves obtained by fitting to experimental data of measured transitions, all calculated by 
Walii et al.l (120151 1. These were combined with a new ab initio dipole moment curve to produce comprehensive line lists for 
these two species. The strong A-X feature about 380 nm should provide a signature for NaH in the atmospheres of cool stars, 
brown dwarfs and exoplanets. However, the shallow nature of the upper A electronic state meansi that this feature is 
predicted to be rather broad. While this means that NaH can make a significant contribution to the opacity, it is likely to 
make this feature difficult to detect. 

NaH was rec ently identified one of only two key diatomic species whose spectral data was completely missing from 
M-dwarf models i Raiourohit et al. 2013tl . A study on the other of these species, AlH, has just been completed and will be 
reported soon. 
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wavelenght, urn 
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Figure 5. Rotational spectrum of ^®NaH computed at T = 298 K; ExoMol (current work) an d CDMS dMu ller et al.ll2005^ . The CDMS 
intensities were obtained by summing in dividual hyperfine t ransitions; these were generated bv lMriller et al.l ll2005tl using an equilibrium 
dipole moment of 6.7 D as calculated bv ISachs et alj (Il975h . 
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